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We study the phase diagram of a symmetric spin-1/2 Heisenberg diamond chain with additional 
cyclic four-spin exchange interactions. The presented analysis supplemented by numerical exact- 
diagonalization results for finite periodic clusters implies a rich phase diagram containing, apart 
from standard magnetic and spin-liquid phases, two different tetramer-dimer phases as well as an 
exotic four-fold degenerate dimerized phase. The characteristics of the established spin phases as 
well as the nature of quantum phase transitions are discussed, as well. 

PACS numbers: 75.10.Pq,75.30.Kz,75.40.Mg 
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I. INTRODUCTION 

Spin systems with cyclic exchange interactions have 
been receiving an increasing amount of attention in 
the past few years. These interactions are known 
to be responsible for the properties of the triangular 
magnetic system composed of 3 He atoms absorbed on 
graphite surfaces. 1 Recently, it was demonstrated that 
a certain amount of four-spin exchange could explain 
the neutron-scattering experiments concerning high-T c 
compounds such as L^CuO^, La6Ca8Cu2404ij£, and 
La4SrioCu2404i4 The latter two are spin-ladder com- 
pounds where charge fluctuations such as cyclic hopping 
processes modify the Heisenberg spin Hamiltonian by 
contributing a four-spin interaction term. In particular, 
it was found that these interactions substantially mod- 
ify the spin triplon gaps and frustrate the formation of 
bound triplon states^ On the theoretical side, in spite 
of the numerous numerical studies predicting rich phase 
diagrams, 5 a number of important questions, concerning 
the type of spin orderings and quantum criticalities real- 
ized by the four-spin exchange, remain unsettled^ 

In this paper we analyze the role of the cyclic four-spin 
exchange interactions in the symmetric diamond chain 
(SDC) which is one of the simplest spin models admit- 
ting multispin cyclic exchange interactions. The model 
is defined by the Hamiltonian (see Fig. [TJ 
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Here Si n , s 2n , and <x n are spin-1/2 operators defined 
on the sites of the n-th elementary cell. In our model, 
the standard cyclic four-spin exchange interactions^ are 
slightly generalized by including its bilinear terms in the 



exchange parameters J, J±, and Jj_. In what follows the 
action and the energy are measured in the units of h and 
Ji, respectively. 

The frustrated diamond chain 8 as well as its various 
modifications, such as the distorted diamond chain^iiS 
and the so-called AB-i ferrimagnetic chain with Ising 
and Heisenberg spins^, have already been discussed in 
the literature in relation to some quasi-lD magnetic 
materials. 12 For the following analysis it is important to 
notice that the cyclic exchange interaction does not vi- 
olate the local symmetry of the Hamiltonian under the 
exchange of the pair of spins (si n , s 2n ) for each diamond 
in the SDC. Thus, in the important case of spin-1/2 op- 
erators Si„ and s 2n , the system is characterized by L 
local good quantum numbers s n = 0, 1 (n = 1, 2, . . . , L) 
related to the composite spins s„ = s ln + s 2 „: s 2 = 
s n {s n + !)■ Using this local symmetry and the stan- 
dard commutation relations for spin operators, the local 
Hamiltonian h n , n +i can be represented in the compact 
form 



h n ,n+l — En + S n • (<7 n + <T n +i) + J n (T n ■ (T n +x 

Tn+l) ■ (2) 



K 



Here e n /J± = s n (s n +l)/2— 3/4 are fixed numbers (—3/4 
or 1/4) for every sector defined as a sequence of the lo- 
cal quantum numbers [si, s 2 , • • ■ , sl], J n = J + K/A — 
s n (s n + l)-K"/2, and {^4, B} is the anticommutator of two 
operators (A and B). 




FIG. 1: The symmetric diamond chain considered in the pa- 
per. The arrows denote the cyclic four-spin exchange interac- 
tion controlled by the parameter K. 

In the parameter space where the ground state is char- 
acterized by s n = (n = 1,2,...,L), the second and 
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fourth terms in Eq. ([2]) vanish and the model is equiva- 
lent to the spin- 1/2 Heisenberg chain with an exchange 
parameter J + K/A. On the other hand, if the ground 
state belongs to the sector s n = 1 (n — 1, 2, . . . , L), Eq. 
© defines the generic model of an alternating mixed-spin 
chain with the simplest three-spin exchange interaction. 
Finally, if the ground state is in the sector [1, 0, 1,0,...], 
the system is reduced to some modification of the well- 
known orthogonal dimer chain^ 



II. PHASE DIAGRAM IN THE PARAMETER 
SPACE (K,J±) 

It is instructive to begin with the phase diagram of a 
single diamond composed of the spins Si, S2, <xi, and 
er 2 . The diamond eigenstates (see Table 1) consist of 
two singlets (Si and S2), three triplets (T\ M , T^, and 
T//; fi = 0,±1), and the quintet (fi = 0,±1,±2). 
In Figure [H the phase boundaries between the single- 
diamond ground states (Si, S2, T%, and T2) are depicted 
by dashed lines. As discussed below, some of these lines 
coincide with the exact phase boundaries of the diamond 
chain, Eq. JT]). 
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TABLE I: Eigenvalues and eigenstates of a single diamond 
composed of the spins si, S2, cri, and a The eigenstates are 
classified according to the following good quantum numbers 
of the single-diamond cluster: the total diamond spin s p , its z 
component Sp, and the diagonal spins s and a [s — si +S2 and 
a = <Ti + o- 2 : s 2 = s(s + 1), cr 2 = a (a + 1)]. The symbols tf* 
((j, — 0, ±) stand for the canonical basic states of the spin-1 
operators, whereas t" denotes the singlet state of two spin- 
1/2 operators. For brevity, only the s z v = components of 
the triplet (T(*) and quintet (Q M ) states are presented. 



A. Phases in the sector [0, 0, 



■ 0] 



For large enough values of the parameter J± , the off- 
chain spins form local dimers, Si n • S2„ = —3/4 so that 
the ground state of the model belongs to the sector 
[0, 0, • • • ,0]. Thus, Eq. © reduces to the form 

(3) 

This is the Hamiltonian of a spin- 1/2 Heisenberg chain 
with the exchange constant J+K/A. The single-diamond 
line AB (defined by J + K/A = 0) coincides with the ex- 
act boundary between the fully polarized ferromagnetic 



phase FM1 ( J + K/A < 0) and the critical spin-fluid 
phase SF ( J + K/A > 0). 
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FIG. 2: Phase diagram of the model, Eq. |TJ, in the param- 
eter space (K, J± ) for J = as obtained from the numerical 
diagonalization of periodic clusters containing up to L = 12 
unit cells. The dashed lines denote the boundaries of the 
single-diamond phase diagram containing the states Si, 1S2, 
Ti, and T2 (see Table 1). The following abbreviations for 
the phases of Eq. |T} are used: critical spin-fluid state (SF), 
fully-polarized ferromagnetic state (FM1), ferrimagnetic state 
(FM2), tetramer-dimer states (TD1,TD2), dimerized singlet 
state (DM), and another ferrimagnetic state (FM3). 

Besides the well-documented collective modes, these 
phases exhibit additional single-particle modes (related 
to the off-diagonal spins) describing excited (si n ,S2 n ) 
dimers, s n = 1. Being eigenstates of the Hamiltonian, 
Eq. ((T|), these excitations are completely localized. Ac- 
cording to Eq. ©, an excited dimer (say at n = L) is 
described by the Hamiltonian 



Hi = C ■ 



K 
~A 



L-l 



C n+1 
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where C 



h,; = 



3 T 
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(L — 1) and hi reads 

J-— I ctl • eri + JtS L ■ ((TL + cri) 



{s L ■ cr L , s L ■ cri}, Ji = 1. 



(5) 



Similar models, describing spin-S* impurities in spin-1/2 
antiferromagnetic Heisenberg chains, have been discussed 
in the literature.— In particular, the above model may 
be considered as a special class of open chains with sym- 
metric couplings of the end spins (<7i and cri/) to the 
external spin (s£ = 1). 

Turning to the region occupied by the SF phase, 
renormalization-group arguments predict that symmet- 
ric perturbations to the open chain are, at most, 
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marginal. Thus, the marginally relevant perturbation 
JiS L ■ (cr L + (Ti) (Jj > 0) is expected to renormalize to 
oo in the SF phase. This corresponds to a fixed point 
where the end spins <Xi and er^ are effectively removed 
from the chain to screen the external spin s^. In terms 
of the original model, Eq. (JTJ, the spins eri, ctl, s± l , 
and S2l form a decoupled single-diamond state (S2) in 
the low-energy sector of the spectrum. Such local excita- 
tions are relevant relatively close to the phase boundary 
AB, where the coupling J^s^ • (a l+ci) dominates the 
biquadratic term in Eq. ([5]). 

For larger values of the parameter K, the biquadratic 
exchange in hi becomes important. Since the energy 
of the quintet state Q grows up with K (see Table 1), 
we concentrate on the triplet state T\ whose energy de- 
creases with K and crosses the energy level of S2 at 
K = 4/5. As discussed in Ref. 14, such a local triplet 
state does not correspond to a stable fixed point, since 
the antiferromagnetic interaction of the effective spin-1 
impurity with the rest of the chain is marginally relevant. 
Thus, one expects that the couplings to the next two 
spins in the chain, <Tl-i and er 2 , grow up to 00. There 
appears another fixed point where the latter two spins 
also decouple from the chain in order to screen the spin of 
the triplet state T\ . Note that in sectors containing more 
spin-1 dimers, other types of screening are possible, as 
well. For example, the spins of neighboring T\ diamonds 
in Fig. [3] may be screened by forming a singlet state. As 
before, the number of decoupled chain spins is four. We 
suggest that such decoupled singlet states, composed of 
longer S2 or Ti diamond chains, control the observed in- 
stability of the SF phase upon decreasing the parameter 
Jj_ . Note that similar complexes of S\ diamonds can not 
produce the instability, since the energy of the resulting 
product state [eri, er 2 ] [er 3 , er 4 ] • • ■ [itl-i, <tl] exceeds the 
energy of the critical phase for arbitrary J±,K > 0. 



B. Tetramer-dimer phases 



close to the AC line separating the single-diamond states 
S\ and 52- On the other hand, the phase boundary be- 
tween the TD1 and FM1 phases clearly deviates from 
the single-diamond boundary AF. The deviation from the 
line AF is related to weak but finite interactions between 
neighboring S2 diamonds appearing in the second-order 
perturbation theory in the parameter J + K/A. 

On increasing the cyclic exchange parameter K at fixed 
J± one finds another tetramer-dimer type state (TD2). 
The picture of the spin-spin correlations, Fig. [31 sug- 
gests that every second diamond is approximately in 
the T\ state. Nevertheless, the numerical analysis shows 
relatively strong antiferromagnetic correlations between 
neighboring T\ diamonds, as opposed to the TD1 state 
where the S2 diamonds are weakly correlated. Clearly, 
both tetramer-dimer phases are gapped and doubly- 
degenerated. According to the general rules^ one may 
expect a discontinuous quantum phase transition be- 
tween the quantum gapped phases TD1 and TD2. Nu- 
merically, the transition point is indistinguishable from 
the exact single-diamond phase boundary CE (K = 4/5) 
separating the single-diamond states £2 and T\. 
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FIG. 3: Short-range spin-spin correlations vs. K in the 
tetramer-dimer phases TD1 and TD2 for periodic chains with 
L = 8 cells (J = 0,J± = 1.15). (1,2) = (cr n ■ s ln ), (1,3) = 
(er n ■ ovfi), (3, 4) = (er„+i • er n+2 ), (2, 5) = (si„ ■ s ln+2 ). 



A detailed numerical study of periodic clusters con- 
taining up to L = 12 cells close to the lines FA, AC, 
and CD in Fig. [2] suggests that for smaller values of J± 
the ground state belongs to the sector [1, 0, • • • ,1,0]. In 
the special case K — 0, our numerical results reproduce 
the phase diagram of the frustrated SDC^ where the so- 
called tetramer-dimer phase (denoted as TD1 in Fig. [2]) 
appears in the interval 0.909 < J± < 2. This doubly 
degenerate singlet state may be roughly represented as 
a product of single-diamond S2 states on every second 
diamond, as depicted in Fig. [3j The short-range corre- 
lations shown in Fig. [3] imply that the simple product 
state is a good variational wave function over the entire 
region occupied by the TD1 phase, excluding a narrow 
region near the line AE. The product state is an exact 
ground state at (K, J±) = (0, 2). An extrapolation of the 
exact-diagonalization (ED) data suggests that the phase 
boundary between the TD1 and SF phases lies extremely 



C. Phases in the sector [1,1, 



■ 1] 



A numerical inspection of the short-range correlators 
(si n • S2„) in finite periodic chains implies that for mod- 
erate values of the parameter J± the established ground 
states belong to the sector [1, 1, • • • ,1]. Thus, in the low- 
energy region the diamond model, Eqs. ([T]) and is 
equivalent to the following mixed-spin Heisenberg model 
(|s n | = Si = 1, |er„| = s 2 = 1/2) with three-spin ex- 
change interactions 



Hi = 51 [s n • (er„ + er n+ i)- J er„ 



K 

1 ~~2<~ Sn ' tr «i S "'' T n + lf 



c n+1 



(6) 
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where J = 3K/4 — J. To the best of our knowl- 
edge, generic mixed-spin quantum Hciscnberg models 
with multiple-spin exchange interactions have not been 
discussed in the literature, although these interactions 
may play an important role in some recently synthesized 
mixed-spin magnetic materials . 16 i 17 Here we restrict our- 
selves to a general overview of the spin phases in the 
specific case (si, S2) — (1, 1/2), and to ferromagnetic ex- 
change interactions between the <x„ spins (J > 0). The 
first two terms in Eq. ([6]) define a standard mixed-spin 
Heisenberg model containing additional non-frustrated 
J exchange bonds. The Lieb-Mattis theorem^ pre- 
dicts a ferrimagnetic ground state for this bipartite model 
which coincides with the classical ferrimagnetic two- 
sublattice Neel state (the FM2 phase in Fig.©. To study 
the role of the competing three-spin interactions, we use 
a qualitative spin- wave analysis supplemented by numer- 
ical ED calculations for finite periodic chains. 



by s„ = si(0,0, 1) and cr n — s 2 (0,0, — 1) survives up 
to k = 1. For k > 1, the ferromagnetic arrangement 
of the smaller S2 spins remains unchanged whereas the 
orientation of the si spins deviates from the z axis: 
s n = si (sin 6 cos <p n , sin 6 sin <fi n) cos 6). Here cos 9 = l/n 
and the azimuthal angle <fi n takes arbitrary values from 
the interval < <j> n < 2ir (n = 1, . . . , L). As a rule, 
quantum fluctuations favor the planar spin configurations 
(<fi n = 0, 7r; n = 1, 2, . . . , L), so that one may expect a re- 
duction of the degeneracy to 2 L . Such high degeneracy 
of the ground state is typical for a number of spin mod- 
els on corner-sharing lattices. As seen in Fig. UJd, in a 
spin-wave approximation the degeneracy produces a full 
line of zero modes lof = in the Brillouin zone. Note 
that the presence of the gapless antiferromagnetic mode 
loaf oc k (k <C 1) for k > 1 is related to the finite trans- 
verse magnetization of the classical canted state. The 
explicit expression for loaf reads 
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FIG. 4: Spin-wave spectrum of the (si,S2) = (1,1/2) chain 
with three-spin exchange interactions, Eq. (|6]), in the ferri- 
magnetic (a) and the canted (b) ground states at k = 0.8 and 
1.5, respectively. J = Jj_ = 0, k = K/2. 
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FIG. 5: Lowest-energy excited states in a periodic chain with 
L = 12 diamonds in the sectors with magnetic moments Mo — 
1 (triangles) and Mq + 1 (circles). Mo — L(Si — S2) = 6 is the 
magnetic moment of the ferrimagnetic ground state. K = 1, 
J = J± = 0. 



Let us start with a discussion of the semi-classical 
limit of the model I©. Since the three-spin interac- 
tion contains an additional overall factor of S1S2, it 
is convenient to redefine the coupling constant K and 
measure the strength of these interactions in terms of 
k = S1S2K. The classical Neel configuration defined 



LOAF = 2J S 2 



£1 

A - 2 



a cos k , (7) 



where a = (k 2 — 1) /(2k J ). It is clear that the antifer- 
romagnetic mode loaf is stabilized by the ferromagnetic 
couplings between the er„ spins. The above picture of 
low-lying excitations in the ferrimagnetic phase close to 
the phase transition point is confirmed by the numerical 
results at K = 1, Fig.0 The numerical estimate for the 
phase transition point is K c w 1.2. 




FIG. 6: Short-range spin-spin correlations vs. K for periodic 
chains with L = 12 cells (J = J± — 0). 



Additional information about the phase diagram of 
Eq. (O may be obtained from the behavior of the short- 
range correlations, Fig. and the singlet-triplet exci- 
tation gap, Fig. with the parameter K. For larger 
K > 1.2, the ED results indicate at least two additional 
phases denoted by DM and FM3 in Fig. [2 DM is a non- 
magnetic gapped singlet state stabilized approximately 
in the region 1.5 < K < 2.3 (see Fig. EJ, and FM3 
is a magnetic phase similar to the FM2 phase. Since 
the unit cell consists of three spin- 1/2 variables, the 
Lieb-Schultz-Mattis theorem 19 suggests that the spin gap 
must be accompanied with at least doubly-degenerate 
ground states. 

Further information about the DM phase can be ex- 
tracted from the dimerization operator D n = S n ■ S n+ i — 
S n • S,i_i where S„ is the spin operator at site n. The 
lattice sites n—1 and n+1 are supposed to be symmetric 
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under the reflection from the central site n. It is conve- 
nient to use the symmetric (antisymmetric) combinations 
\s, a) — (|0) ± |l))/\/2, where |0) is the translationally in- 
variant singlet ground state of the finite periodic chain 
and 1 1) is the singlet excited state which is almost degen- 
erate with |0). The states \s,a) are not translationally 
invariant and may be expected to produce finite values of 
the dimer order-parameter (s,a\D n \s,a) in a dimerized 
system.— The caricature of the DM state presented in 
Fig. [5] is obtained from the extrapolation of the ED re- 
sults for (s, a\D n \s, a) (L = 8, 10, and 12). A pronounced 
enhancement of the dimer order-parameter with L is in- 
dicated only for two types of bonds in the dimer model 
(the thick lines in Fig. [5]) . Clearly, the suggested ground 
state exhibits a four-fold degeneracy. 
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FIG. 7: The excitation gap in the DM phase as a function of 
the parameter K. 

Finally, the numerical results point towards the exis- 
tence of another partially polarized magnetic phase in the 
narrow interval 1.2 < K < 1.5 between the ferrimagnetic 
(FM2) and the dimerized (DM) phases (the hatched area 
in Fig. [5]). Recently, similar exotic magnetic states have 
been predicted in a number of one-dimensional spin sys- 
tems with magnetic frustrations Q Typically, discussed 
spin states exhibit a partially polarized magnetization in 
the z direction (M < M = s\ — s 2 ), a quasi- long-range 



transverse magnetic order, and a gapless linear mode re- 
lated to the destroyed classical canted state. Unfortu- 
nately, the methods used in the present study do not 
suggest a clear statement indicating the existence of such 
an exotic phase in the discussed system. 



III. SUMMARY 

In conclusion, we have examined the impact of the 
cyclic four-spin exchange interactions on the ground- 
state phase diagram of the symmetric spin-1/2 diamond 
chain. Using the local symmetries of the model, the spin 




FIG. 8: Caricature of the dimerized singlet phase DM as ob- 
tained from an extrapolation of the numerical ED results for 
periodic chains with L = 8, 10, and 12 . The open circles 
(thick lines) denote the symmetric (antisymmetric) eigen- 
states of two spin-1/2 operators. K — 2, J± — J = 0, 

(8 U 82) = (1,1/2). 



phases were classified by the set of good quantum num- 
bers s„ = 0, 1 (n — 1, 2, . . . , L) related to the composite 
spins s„ = Si„ + S2„. The presented study demonstrates 
a rich phase diagram in the parameter space (K,J±). 
Apart from the standard magnetic and paramagnetic 
phases, the system exhibits two different tetramer-dimcr 
phases in the sector [1, 0, 1, 0, ■ • • ,1,0] as well as an ex- 
otic four-fold-degenerate dimerized ground state in the 
sector [1, 1, • • • ,1]. 
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